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Abstract 

A new type of deterministic chaos for a system described by iterative two-dimensional maps is 
reported. The series being generated by the original map has an average upward trend while the 
first difference, which is the series of changes from one period to the next, exhibits chaotic behavior 
through period doubling bifurcation. Furthermore, step-wise time series appears as the values of 
the bifurcation parameter are large, and the first difference of the time series exhibits intermittent 
chaos. 
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In 1974 Robert May [1,2] first reported on a one-dimensional map, such as the logistic 
map, which generates a chaotic dynamics. He proposed the non-linear difference equations as 
a model of population growth which takes place at discrete intervals of time and generations 
are completely non-overlapping. The time series appeared from the low-dimensional maps 
proposed by May is stationary. Namely, its statistical properties such as mean, variance, etc. 
are all constant over time. Therefore, theorems on existence of chaotic attractors generated 
by iterative maps [3,4] are necessarily based on the assumption that the time series can be 
rendered stationary. Since their studies, it has been believed that re-injection mechanisms 
to suppress exponential divergence of the variable are indispensable to emergence of chaos. 

The aim of this letter is to show a new type of deterministic chaos for an iterative two- 
dimensional map which generates non-stationary time series. More concretely, the time series 
x(t) being appeared from the original map which we here propose has an average upward 
trend while the first difference of the time series of x(t), which is the series of changes from 
one period to the next, Ax(t) = x(t) — x(t — 1) exhibits chaotic behavior through period 
doubling bifurcation. Furthermore, as the values of the bifurcation parameter are large, 
step-wise time series of x(t) appear and the first difference of x(t) exhibits intermittent 
chaos. 

Consider the two-dimensional difference equations 

x(t + 1) - x(t) = a[exp(y(t) - x(t)) - 1], (1) 

y(t + l)-y(t)=0[x(t)-y(t)]. (2) 

Equation (2) can be transformed into an alternative form 

oo 

y (t) = / 3^(i-/?r i x(t-2). (3) 

1=1 

The variable y{t) is an exponentially weighted moving average of x(t). The weight (1 — /?)* _1 
is assumed to be exponentially decreasing, that is, < j3 < 1. After this, we specify the 
parameters as follows : a > 0, and (3 = 0.5. Under these conditions, x(t) tends to go back 
to the exponentially weighted moving average, y(t). Namely, if y(t) > x(t), then x{t + 1) 
increases, and y(t + 1) decreases. The dynamic system (l)-(2) has infinite number of fixed 
points: x(t) = y(t). As can be checked easily, the local stability condition of a fixed point is 
a + /3 < 2. The time series of x(t) has an average upward trend for the condition, a + /3 > 2. 



The upward trend is robust with respect to the initial values of x(t) and y(t). Figure 1 is 
a sample path of x(t) as a function of time for a = 5 and /? = 0.5. The time series of x(t) 
fluctuates around the upward trend. To stationarize a time series of x(t), we take the first 
difference Ax(t) of a time series of x(t) which is the series of changes from one period to the 
next. Figure 2 shows that the corresponding time series of Ax(t) for a = 5 and (3 = 0.5. The 
sample path of the first difference Ax(t) in Figure 2 is stationary and chaotic. To account 
for the chaotic behavior, let us consider the first return map of the form: 

Ax(t+l) = f(Ax(t)). (4) 

Figure 3 is the bifurcation diagram for the first difference, Ax(t) generated by the first 
return map f(Ax(t)) where the parameter a varies smoothly from 1 to 5. Figure 3 shows 
that chaos occurs through a sequence of period doubling as the parameter a increases. To 
construct the first return map, Figure 4 plots Ax(t) as a function of Ax(t) which is obtained 
for a chaotic regime of the attractor with a = 5 and f3 = 0.5. One can reconstruct the graph 
of the first return map f(Ax(t)) from these observations. Figure 4 shows that we obtain a 
curve which has an extremum. 

The bifurcation diagram shows the period-three attractor located at a = 6. Figure 5 
show that the first retuen map f(Ax(t)) has a period-three point for a ~ 2.9 and j3 = 0.5 
The Li- York theorem [3] implies that such a map has points with every possible period. 

As a increases further, step-wise time series of x(t) appears. Figure 6 is the time series 
of x(t) for a = 15. Figure 7 is the first difference Ax(t) of the time series of x(t) in Figure 
6. The time series of Ax(t) exhibits intermittency. As the parameter a is large, the laminar 
phase duration of intermittency is longer. 

In this letter we propose a dynamical model which exhibits a new type of non-stationary 
chaos. Our model can be widely applied to non-stationary time series having an upward 
trend such as human population, World GDP, and so on. 
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FIG. 2: The time series of Ax(t) = x(t + 1) — x(t) as a function of time t for a = 4, and /3 = 0.5. 
The time series corresponds to the first difference of the time series of x(t) in Figure 1. 
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FIG. 3: The bifurcation diagram for the first difference of the time series of x(t), where the 
parameter a varies smoothly from 1 to 5. 
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FIG. 5: The period-three of the first return map for a = 2.9 and (3 = 0.5. 
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FIG. 7: Intermittency: The time series of Ax(t) = x(t + 1) — x(t) for a = 15, and (3 = 0.5. The 
time series corresponds to the first difference of the time series of x(t) in Figure 6. 



